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ABSTRACT 

The rotation curves and the relative mass distributions of the two nearby Local 
Group spiral galaxies, M31 and M33, show discrepancies with Modified Newtonian 
dynamic (MOND) predictions. In M33 the discrepancy lies in the kinematics of the 
outermost regions. It can be alleviated by adopting tilted ring models compatible 
with the 21-cm datacube but different from the one that best fits the data. In M31 
MOND fails to fit the falling part of the rotation curve at intermediate radii, before 
the curve flattens out in the outermost regions. Newtonian dynamics in a framework of 
a stellar disc embedded in a dark halo can explain the complex rotation curve profiles 
of these two galaxies, while MOND has some difficulties. However, given the present 
uncertainties in the kinematics of these nearby galaxies, we cannot address the success 
or failure of MOND theory in a definite way. More sensitive and extended observations 
around the critical regions, suggested by MOND fits discussed in this paper, may lead 
to a definite conclusion. 

Key words: galaxies: halos galaxies: individual: M31,M33 galaxies: kinematics and 
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1 INTRODUCTION 

According to Newtonian dynamics the mass distribution of 
the luminous components of spiral galaxies cannot account 
for the observed profiles of their rotation curves (hereafter 
RC) especially in the outer regions of galaxies (Rubin, Ford, 
& Thonnard 1980). To account for this discrepancy, which 
becomes more pronounced for late type galaxies, discs are 
thought to be embedded in dark halos of non baryonic mat- 
ter (Persic, Salucci & Stel 1996). An alternative explanation 
for the mass discrepancy has been proposed by Milgrom 
by means of the modified Newtonian dynamics or MOND 
(Milgrom 1983). According to this theory the dynamics be- 
comes non-Newtonian below a limiting acceleration value, 
ao ~ 10 -8 cm s -2 , where the effective gravitational acceler- 
ation takes the value a e // = yJaog n , with g n the acceleration 
in Newtonian dynamics. Outside the bulk of the mass dis- 
tribution, MOND predicts a much slower decrease of the 
(effective) gravitational potential, with respect to the New- 
tonian case. This is often sufficient to explain the observed 
non-keplerian behavior of RC (Sanders & McGaugh 2002, 
Lokas 2002). 

This success is remarkable in that MOND has only one 
free parameter, namely ao. MOND theory of Gravity has 
been generalized in the framework of a relativistic, Lorentz 
invariant theory (Bekenstein 2004). It is worth noticing that 
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MOND could account for the acoustic peak structure of Cos- 
mic Microwave Background experiments (Slosar, Melchiorri, 
& Silk 2005, McGaugh 2004). Gravitational lensing data to- 
gether with X-ray data in galaxy cluster regions give perhaps 
the most convincing proof for the existence of "dark matter" 
(Clowe, Gonzalez & Markevitch 2004; Pointecouteau & Silk 
2005) though on much larger scales than we consider here. 
On the other side, there are difficulties for the dark matter 
paradigm as well (e.g. Goerdt et al. 2006) and dark particles 
have not been detected yet. 

It is however important to check MOND validity on 
galactic scales, that host the empirical phenomenon which 
stimulated its birth and its theoretical development. MOND 
is unable to fit some RC of spiral galaxies (e.g. Gentile et al. 
2004) but this failure is not decisive to totally rule out this 
theory of Gravity. In fact MOND mass modeling depends 
on the exact value of the acceleration threshold ao, on the 
galaxy distance and on the detailed distribution of stars and 
gas. Uncertainties in these quantities weaken the relevance 
of the comparison between model predictions and data. All 
the above, in addition to the possible presence of bars, in- 
teractions and warps, makes the MOND-dark matter debate 
still open (Bottema et al. 2002). In this perspective M31 and 
M33, well studied spirals and primary distance indicators, 
are excellent benchmarks for MOND. For these galaxies, all 
the above uncertainties are fairly small: distances are well 
known (780 and 840 kpc within 10%), and their exponential 
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stellar discs and gas surface brightness are precisely mea- 
sured. 

M33 is one of the fewest objects for which it is possible 
to combine a high quality and high resolution RC from CO 
J=l-0 line that extends inward to 200 pc with a high qual- 
ity 21-cm line RC which extends out to 13 disc scalelengths 
i.e. 19 kpc (Corbelli & Salucci 2000, Corbelli 2003). The 
galaxy has no prominent bar or bulge and the HI and CO 
velocity fields are very regular and cannot be explained by 
Newtonian dynamic without including a massive dark mat- 
ter halo. M31 is a very large nearby galaxy whose extended 
and spatially resolved RC has been recently studied, in the 
Newtonian framework, to advocate a dark matter compo- 
nent contributing more than half the mass inside 30 kpc 
(Widrow, Perrett & Suyu 2003, Geehan et al. 2006, Carig- 
nan et al. 2006). 

In this paper we use recent kinematical data of these 
two galaxies to test MOND theory. It is important to in- 
vestigate the uncertainties left in the kinematics or mass 
distribution in order to claim complete failure or success of 
this theory. To this purpose we investigate M33 in Section 2, 
and in Section 3 we check the compatibility of MOND with 
less probable but still possible deconvolution models for the 
M33 21-cm datacube. In Section 4 we analyze MOND mass 
models for M31 using the most recent RC data. Section 5 
summarizes our results and addresses future observational 
strategies. 



2 THE M33 ROTATION CURVE IN THE 
MOND FRAMEWORK 

Neutral gas in M33 extends well beyond the optical disc. The 
rotation curve and the average radial distribution of gas can 
be derived from 21-cm spectral line data using a set of tilted 
concentric rings that accounts for the possible presence of 
a warp. This is especially important when testing MOND 
because the kinematics and the distribution of matter in 
the outermost regions of galaxies are strongly coupled by 
the local non-Newtonian gravitational field. 

The RC that we use is shown in Figure 1 with 2a er- 
rorbars in each bin derived from the dispersion of single 
spectra velocities after deconvolution. Data have been de- 
convolved using the tilted ring model which best fits the 
Arecibo 21-cm data (hereafter deconvolution Model 1, see 
bottom panels of Figure 1; Corbelli & Schneider 1997, Cor- 
belli & Salucci 2000). Uncertainties in this model will be 
investigated in the next Section. The inner RC, over the 
interval 0.2 < R < 3.4 kpc, is from CO J=l-0 data at 
an angular resolution of 0.75 arcmin (Corbelli 2003). For 
3.4 ^ R < 5.5 kpc we use 21-cm interferometric observa- 
tions of Newton (1980) at 1.5x3 arcmin angular resolution. 
The outer RC is from Arecibo 21-cm data at 4.5 arcmin 
angular resolution (Corbelli & Schneider 1997). Corbelli & 
Salucci (2000) and Corbelli (2003) have used these data for 
modeling the M33 mass distribution when a dark matter 
halo is in place. Note however that the distance assumed in 
Corbelli (2003) is 760 kpc and not 840 kpc as erroneously 
stated in that paper. Here we shall use a distance to M33 
of 840 kpc (Freeman et al. 2001). As shown by Corbelli & 
Salucci (2000) the rotation curves of the receding and ap- 
proaching sides of M33 look very similar for deconvolution 



Model 1. Their slopes are consistent (within la errors) with 
those relative to the rotation curve of the galaxy as a whole. 

The visible mass components in M33 are: stars, dis- 
tributed in a disc and in a spheroid, and gas, in molecular 
and in neutral atomic form. In the stellar disc we assume 
that the mass follows the light distribution in the K-band 
that is well fitted by an exponential law with scalelength 
R d = 5.8±0.4 arcmin (1.4±0.1 kpc) (Regan & Vogel 1994). 
There are no wiggles or deviations from a pure exponential 
disc between 3 and 18 arcmin. Beyond 18 arcmin the the 
near-infrared brightness drops below the sensitivity limit of 
the quoted observations. However preliminary work in the 
I-band by Ferguson et al. (2006) shows that the exponen- 
tial surface brightness holds out to about 40 arcmin. In this 
Section we extrapolate such exponential stellar light distri- 
bution beyond 40 arcmin, to the outermost RC data point. 
Alternative models are discussed in Section 3. For consis- 
tency with other works on MOND we assume an infinitely 
thin stellar disc. In the innermost region (R < 3 arcmin) 
there is a luminosity excess with respect to the exponen- 
tial disc extrapolation and a small spheroid is likely to be in 
place though its luminosity and mass are very uncertain (Re- 
gan & Vogel 1994; Corbelli 2003). We will take this stellar 
component into account by parameterizing its contribution 
to the circular velocity as: 

V sph = GM aph R"- 5 /(R + s) (1) 

the total mass of the spheroid, M ap h, and its typical ra- 
dius, s, are free parameters. Observations suggest values of 
s in the 0.1-1 kpc range and spheroidal masses smaller than 
10% of the stellar disc mass (see Corbelli 2003 and refer- 
ences therein). Being M33 a blue galaxy, also in its central 
region, the only restriction we place on the spheroid here 
is M B ph/L S ph < 3 where L sp h = 4 x 10 8 L is the largest 
blue spheroidal luminosity ever estimated. The atomic and 
molecular gas surface densities in M33 are shown by Corbelli 
(2003). They result from azimuthal averages of the data af- 
ter deconvolution according to the tilted ring model which 
best fits the 21-cm datacube. The molecular mass of M33 is 
less than 10% the HI mass using the CO to H2 conversion 
factor determined by Wilson (1995) in M33. 

The free parameters of the M33 mass model in the 
MOND framework are: s, M ap h, and Md/Ld, (where Ld = 
5.7 x 10 9 Lq is the blue disc luminosity). We use the critical 
acceleration value ao derived from the analysis of a sam- 
ple of rotation curves ao = 1.2 x 10~ 8 cm s~ 2 (Sanders & 
McGaugh 2002). 

We use the reduced chi-square statistic, xt, to judge the 
goodness of a model fit. Figure 1 shows the MOND best-fit 
model curve: its xi = 2.0 is higher than the value 0.7 ob- 
tained for dark halo models (Corbelli 2003). The best fitting 
values of the free parameters are Md/Ld = 0.6, s = 1.4 kpc, 
and M aph = 1.2 x 10 9 M (the maximum allowed value). 
The fit is poor and only marginally acceptable, being 99.9% 
the probability of finding \v ^ 2.2 for random events. The 
rotational velocities predicted by MOND are higher than ob- 
served for R ~ 7— 8 kpc and lower than data for R > 15 kpc. 
MOND has difficulties to reproduce the M33 RC because of 
the rising rotational velocities in the outer disc where the HI 
gas mass density strongly declines. Moreover, the resulting 
spheroidal mass, about 20% of the total stellar mass, is un- 
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Figure 1. Data for M33 rotation curve from deconvolution Model 
1 and the MOND best fit. The long and short dashed lines are 
the Newtonian RC of the spheroid and stellar disk respectively. 
The dot-dashed line is the gas Newtonian RC. The bottom panels 
show the inclination i and position angle 9 of Model 1. 



realistic for this blue galaxy, it bounces to the extreme value 
allowed by our fit and its typical radius is larger than obser- 
vational estimates. Variations of the stellar disc scalelength 
within the uncertainties quoted by Regan & Vogel (1994) do 
not affect the above conclusions significantly. 

In modeling the gas distribution in the outer disc we 
have assumed that most of the gas is in neutral form and is 
traced by the 21-cm line. Absorption measurements of the 
21-cm radiation from background sources behind the M33 
disc have excluded the possibility that there are relevant 
HI masses in the outer disc of M33 which are undetected 
in emission at 21-cm (Corbelli & Salpeter 1993a). There is 
however the possibility that extragalactic UV radiation ion- 
izes the outer disc gas as the HI column density approaches 
the value of 10 19 cm -2 (Corbelli & Salpeter 1993b). This 
implies a sharp HI radial fall-off, as that observed along the 
M33 major axis by Corbelli, Schneider & Salpeter (1989). 
The gas extends beyond 18 kpc but it is undetected via 21- 
cm radiation because it is fully ionized. The sharpness of the 
HI edge is a function of the gravity perpendicular to the disc 
and it is worth noticing that MOND models reproduce the 
observations as the Newtonian dynamic does in the pres- 
ence of a dark halo. The presence of an ionized component 
extending radially further out would rise the gas RC by 1-2 
km s _1 at the outermost sampled radius. This is however 
not sufficient to improve the significance of MOND fit to the 
RC. 



3 CAN OBSERVATIONS AND MOND 

PREDICTIONS IN M33 REACH A BETTER 
AGREEMENT? 

In this Section we investigate whether MOND models can 
provide a better fit to the M33 RC when uncertainties in 
the modeling presented in Section 2 are taken into account. 

(a) Distance uncertainties 

For M33 we are using a distance of D = 840 kpc de- 
termined from Hubble Space Telescope Cepheid data (Free- 
man et al. 2001). The uncertainties on this measured value 
are ±40 kpc. To be conservative we take into account other 
values obtained by different distance indicators. These fall 
in the range 790-940 kpc (McConnachie et al. 2004; Mc- 
Connachie et al. 2005 and reference therein; Bonanos et al. 
2006; Sarajedini et al. 2006). We model fit the M33 RC by 
assuming the distance as an additional free parameter. The 
best fit is obtained for a distance to M33 of 860 kpc, but the 
improvement in term of the xt value is negligible (1%). 

(b) Departures of the stellar mass distribution from a simple 
exponential law 

Preliminary work by Ferguson et al. (2006) shows that 
the surface brightness of the M33 disc has no departures 
from a pure exponential law out to about 40 arcmin. Be- 
tween 40 and 60 arcmin, only RGB star counts are available 
and their distribution can be fitted by a steeper exponen- 
tial law than that describing the K and I-band inner sur- 
face brightness. If the distribution of RGB stars reflects the 
stellar mass distribution in the outer disc we have to modify 
the stellar mass exponential scalelength used in the previous 
Section at large radii. We shall use a scalelength of 1.4 kpc 
(Regal & Vogel 1993) for R < 10 kpc and of 0.9 kpc for 
10 < R < 20 kpc. For sake of completeness we will consider 
also models with a larger stellar scalelength, up to 2 kpc, for 
R > 10 kpc. 

In any case (for both smaller and larger scalelength at 
large radii) variations of the best fitting parameters and xi 
are hardly noticeable. This is because the change in scale- 
length takes place in a region where the stellar contribution 
to circular velocity is not prominent and radially declining. 

(c) Tilted ring models 

The modest inclination of M33 with respect to our line 
of sight implies that small variations of the inclination angle 
give non negligible variations of the deprojected rotational 
velocity. In this Section we check if tilted ring models, with 
less statistical relevance but still within the 95% confidence 
level for fitting the 21-cm Arecibo datacube, yield a rotation 
curve in closer agreement with MOND predictions. 

The tilted ring model we use has 11 free rings. The 
fitting program interpolates between the 11 free rings to 
find the parameters at each radius (see Corbelli & Schnei- 
der 1997). Here we vary the parameters of the tilted ring 
model as follows. We fix the ring parameters to the best fit- 
ting values of Model 1 for R < 6 kpc, and in the interval 
6 ^ R ^ 15 kpc we vary the inclinations of the 5 free rings 
by ±5deg respect to Model 1 values. Each of these models 
is compatible with the 21-cm datacube and yields a rota- 
tion curve over 0.2 ^ R ^ 15 kpc that we fit using MOND 
prescription. 

The use of deconvolution models whose ring inclina- 
tions are systematically lower than those of Model 1 for 
10 ^ R ^ 15 kpc, results in a RC which is flatter than that 
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Figure 2. Data for M33 rotation curve from deconvolution Model 
2 and the MOND best fit to it. The long and short dashed lines are 
the Newtonian RC of the spheroid and stellar disk respectively. 
The dot-dashed line is the gas Newtonian RC. The bottom panels 
show the inclination i and position angle 8 relative to Model 2. 
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Figure 3. Unbinned kinematical data from 21-cm line observed 
at Arecibo using deconvolution Model 2. Filled triangles refer 
to the approaching side of the galaxy, and open squares to the 
receding side. Straight lines are the best linear fit for R > 6 kpc 
to all data (continuous line), to the approaching side set (dashed 
line) and to the receding side set (dashed-dotted line). 



of Figure 1 and in better agreement with MOND model. We 
then consider for these 5 free rings the deconvolution model 
which gives the lowest xl if° r MOND fit out to R = 15 kpc, 
and we vary the position angle and inclination of the two 
outermost free rings by ±5, 10, 15 deg. We choose a combi- 
nation of position angles and inclinations of the last two 
rings, compatible with 21-cm data, whose associated RC 
has the lowest xl for MOND model fit ( X l = 0.84). Fig- 
ure 3 shows the whole M33 rotation curve derived according 
to the above prescriptions and the relative deconvolution 
model (hereafter deconvolution Model 2). 

MOND fit improves considerably when outermost orbits 
are more inclined and at more negative position angles than 
in Model 1. The velocity fields of the two separate halves of 
M33 look very similar also for Model 2 and are shown in Fig- 
ure 3. The linear fit to the data for R > 6 kpc have consistent 
slopes (within 1-cr error). The parameters of the mass model 
are M d /L d = 0.7, s = 0.6 kpc and M sph = 4 x 10 s M©. 
The MOND best fit to Model 2 data is good and requires 
spheroidal parameters closer to values suggested by multi- 
wavelengths observations of M33. Deconvolution Model 2, or 
similar ones, are therefore acceptable alternatives to Model 1 
and give rotation curves compatible with MOND predictions 
with high statistical significance, similar to dark matter halo 
models. Only deeper 21-cm maps or Ha emission searches 
around and beyond the 18 kpc radius can constrain the de- 
convolution model further and solve the ambiguity left in 
M33 by the current data. 



4 THE M31 ROTATION CURVE IN THE 
MOND FRAMEWORK 

Recent 21-cm observations of M31 along its major axis have 
shown that the rotation curve of this nearby galaxy stays 
flat out to 35 kpc (Carignan et al. 2006). Because of the 
discrepancies between different determinations of the RC in 
the inner 10 kpc and the relative unimportance of this region 
for MOND tests, fully into the Newtonian regime, we con- 
sider only the radial range and HI rotational velocities given 
by Table 1 of Carignan et al. (2006). For this radial range 
(10-35 kpc), there are very few uncertainties on the galaxy 
inclination and position angle. The warp of M31 is a minor 
one, due perhaps to the more massive and extended stellar 
disc. Models of the warp at large radii imply a slightly in- 
creasing inclination as a function of radius (Newton & Emer- 
son 1977; Briggs 1990). Due to the nearly edge on position 
of the optical disc, the use of a constant inclination of 77° 
instead of a higher one at large R overestimate the rota- 
tional velocities by only 2-3%. 21-cm maps suggest higher 
inclinations at large R (Brinks & Burton 1984) but analysis 
of the extended stellar disc by Ibata et al. (2005) suggest 
instead a somewhat lower values of i. The large errorbars 
in the rotation curve of this region take these uncertainties 
into account. 

We consider a mass model with three components: a 
gaseous disc, a stellar disc and a central spheroidal or stel- 
lar bulge. This last component is very prominent and its 
mass can be of the same order of the disc mass. Walterbos 
& Kennicutt (1987,1988) derive the bulge and disc optical 
luminosity profiles. There are very large uncertainties on the 
bulge scalelength and on its mass because the disc-bulge de- 
composition is not very robust. The bulge scalelength range 



Testing MOND with Local Group spiral galaxies 5 



300 
200 

E 

> 100 



300 
200 

6 

> 100 



r 



(a) 



H — I — I — I — KH- 



H 1 h 



t \ N 



(b) 



10 20 
R(kpc) 



30 



Figure 4. Data for M31 rotation curve and the MOND best fit 
to it. The long and short dashed lines are the Newtonian rotation 
curves of the bulge and stellar disc respectively. The dot-dashed 
line is the gas Newtonian rotation curves. Panel (a) shows the 
best MOND fit for r b = 0.61 kpc, panel (6) for r b = 1.8 kpc. 



is 0.61 ^ r b ^ 1.8 kpc (Geehan et al. 2006) and here we shall 
consider the two extremes of this range. For the disc scale- 
length we shall use that K-band photometric data which 
give a value of 4.5 kpc (Battaner et al. 1986). Similar values 
are found in the more extended R-band photometric maps 
(5.1 kpc out to 40 kpc, Ibata et al. 2005). The bulge and 
disc blue luminosities are 9 x 10 9 L Q and 2 x 10 10 L Q re- 
spectively. The observed colors and population models imply 
stellar mass to light ratios between 2.8 and 6.5 in solar units 
(Bell & de Jong 2001) which we will use to limit our models. 

We adopt the neutral gas surface density measured by 
Sofue & Kato (1981). The molecular gas mass is less than 
10% of the gas mass (Nieten et al. 2006) and its peak is 
located at about 11 kpc. Results are insensitive to the inclu- 
sion of this components as well as to the central supermas- 
sive black hole whose estimated mass is ~ 0.5 - 1 x 10 s M 
(Salow & Statler 2004; Bender et al. 2005). 

For r b = 1.8 kpc the minimum xt is 10.3 with Md/Ld = 
4.5 and Mi/Lb = 6.5. For r b — 0.61 kpc the minimum xt 
is 8.1 with Md/Ld = 3.9 and M b /L b = 6.3. MOND fits are 
shown in Figure 4. The evident failure is due to the declining 
part of the rotation curve of M31 between 10 and 20 kpc. 
While the classical Newtonian dynamic with a dark matter 
halo is able to provide a good fit to the data (Widrow et 
al. 2003; Carignan et al. 2006) by making the stellar disc 
responsible for the declining part of the curve and the dark 
halo responsible for the outermost flattening, MOND is un- 
able to reproduce this non monotonic behavior. It is in fact 
around 10 kpc that g n ~ ao and non-Newtonian corrections 
start to be important and force a falling Newtonian RC into 
a flat one, inconsistent with the data. 

Uncertainties on the M31 distance are much smaller 
than for M33. Distances resulting from different data sets 



and methods are always in the range 770-795 kpc. There- 
fore possible variations with respect to the assumed one 
(D=780 kpc) are less than 2% and give negligible effects 
on the MOND fit (Freeman et al. 2001; Joshi et al. 2003; 
McConnachie et al. 2005 and reference therein; Ribas et al. 
2005). Concerning the uncertainties on the disc scalelength 
we have found no evidence for a scalelength smaller than 
indicated by K-band photometry. However from measure- 
ments of light profiles in other bands, extending at larger 
radii, we cannot exclude a flattening of the stellar light pro- 
file at large radii (Walterbos & Kennicutt 1988; Ibata et al. 
2005; Irwin et al. 2005). MOND models which take into ac- 
count a change in the disc scalelength at large radii do not 
give acceptable fits to M31 RC data. 



5 CONCLUSIONS 

Local Group galaxies benefit from low observational uncer- 
tainties in distance determination and in the distribution 
of baryonic matter. Moreover, they provide RC data to an 
unbeatable spatial resolution. For this reason we have cho- 
sen the two nearest spiral galaxies, M31 and M33, to test 
MOND. M33, a low luminosity late galaxy, has a very ex- 
tended HI disc providing a rotation curve out to 13 disc 
scalelengths. Dark matter dominates over the visible bary- 
onic matter well inside its optical radius according to New- 
tonian dynamic. However its inclination is not sufficiently 
high to limit the uncertainties on the rotational velocities 
of the outer disc. Here the presence of a warp leaves some 
degree of freedom to the tilted rings deconvolution model. 
Even though MOND models are marginally compatible with 
the rotation curve from the best fitting tilted ring model to 
the 21-cm datacube, data leave some space for a less ris- 
ing curve, well fitted by MOND models. Since deconvolu- 
tion models in better agreement with MOND theory predict 
a more inclined outermost disc, future observations should 
check whether this is effectively the case by searching for 
faint emission above the main disc at a specific velocity. Sen- 
sitive searches for extended emission should be carried out 
both at 21-cm and in the Ha line (since there is the possi- 
bility that extragalactic UV radiation ionizes the outermost 
disc) . 

In M31 the disagreement between the RC and MOND 
predictions is statistically significant. MOND models are un- 
able to reproduce the non monotonic behavior of the rota- 
tion curve in regions where g„ ~ ao and non-Newtonian cor- 
rections predict a flat RC. However the actual RC is made up 
of data with small uncertainties in some radial range while 
in others and beyond 20 kpc uncertainties are large. For this 
galaxy it would be nice to have more uniform uncertainties 
by producing a fully sensitive two dimensional map at 21-cm 
which provides kinematical information on a more extended 
area and defines the orientation of the possible warp. 

We are grateful to the anonymous referee for her/his 
useful criticism to the original manuscript which have im- 
proved the quality of the paper. 
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